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图 论 中 最 短路 问题 的 MATLAB 程序 实现 
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搞 要 :; 解决 图 论 中 最 短路 问题 的 最 好 方法 --“Dijstra 算法 ”, 通 过 解析 实例 模型 ,对 模型 算法 进行 描述 .拓展 ,并 给 出 
了 求 最 短路 以 及 求 最 短路 长 的 MATLAB 程序 ,此 程序 具有 通用 性 。 
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1 问题 的 提出 
设 CE = (Vy, 局 为 连通 图 ,顶点 集 为 {1, 2, 3,…nj, 图 中 各 边 (i, ) 有 非 负 权 ci ( 当 (i,j) 不 是 边 时 , 权 等 
于 inf; 当 Gj) 是 有 向 边 时 ,ci 与 0; 可 以 不 相等 ), 求 一 条 道路 使 它 是 从 顶点 1 到 顶点 n 的 所 有 道路 中 总 
权 数 最 小 的 路 ,这 就 是 图 论 中 的 最 短路 问题 mn。 解 决 这 个 问题 至 今 公认 最 好 的 方法 是 1959 年 提出 Di- 
jkstra 算法 , 它 用 于 计算 一 个 点 * 到 其 他 所 有 点 的 最 短路 。 此 算法 基本 原理 是 : 若 某 条 路 是 最 短路 ,这 
条 路 上 的 任意 一 段 路 也 是 连接 这 段 路 两 个 端点 的 最 短路 "3。 
2 算法 描述 84 
第 一 步 将 点 标 上 永久 性 标记 Ps)=0, 表示 从 开始 点 s 到 达 上 点 s 的 最 短 距 高 是 零 。 
第 二 步 将 其 余 的 顶点 标 上 了 标记 ,T 记号 是 试探 性 标记 ,点 j 的 T 标 记 70) 分 两 部 分 ,TO)=T,()(T0))， 
第 一 部 分 TO) 为 从 开始 点 经 过 带 P 标 记 的 点 到 达 j 点 的 最 短路 的 权 和 ,括号 中 720) 为 第 二 部 分 ,是 这 
最 短路 中 j 点 的 前 一 点 (如 有 多 条 最 短路 , 则 720) 可 能 有 多 值 )。 不 能 经 过 带 P 标 记 的 点 到 达 的 点 的 工 值 
是 无 穷 大 (用 inf 表示 ), TT 是 空 集 。 
第 三 步 若 这 些 带 有 了 标记 中 权 和 数 7 最 小 的 点 是 , 则 点 是 带 忆 标记 的 点 外 与 开始 点 s 最 近 的 - 
点 。 把 点 有 的 TT 标记 改 为 PP 标记, 如 果 权 和 数 最 小 的 点 有 多 个 , 则 把 它们 都 改 为 P 标 记 。 若 点 nn 不 是 P 
标记 , 转 第 二 步 ( 对 带 有 7 标记 的 点 重新 标记 ,直至 点 4 为 标记 为 止 )。 
第 四 步 追寻 最 短路 , 从 终点 n 开始 逆向 逐次 求 最 短路 经 过 的 点 权 和 为 Pn). 
从 算法 直接 可 见 所 得 到 的 路 是 最 短路 。 上 述 算法 更 具体 的 步骤 如 下 :不 妨 设 开始 点 的 标 导 是 1。 
() 设 N = 0，Pl) = 0, 其 余 各 点 都 是 了 标记 ,mh 值 为 无 穷 大 , 7 值 为 空 集 。 
(2) 若 w 点 为 刚 成 为 已 标记 的 (一 个 或 几 个 ) 点 ,将 所 有 与 这 些 w 相 邻 的 带 有 了 标记 的 点 v 的 了 标 
记 的 值 改 为 Ti(Yj,N +1)= min [TY), N), Py)+ Cs); 
车 TilosN+l) = Po + cp 则 To， N+1) =v; 车 TiopN+l = TsN, 则 Tvj ,N+1) =7ox N，( 因 此 
ZT 可 能 是 多 值 的 )。 | 
(3) 比 较 所 有 具有 了 了 标记 的 点 ,把 其 中 思 值 最 小 者 的 点 内 的 了 标记 改 为 已 标记 ,Pon = minTi(vj, N+ 
1); 当 存 在 两 个 以 上 的 最 小 者 时 ,可 同时 将 它们 改 为 P 标 记 ; 车 点 w = n 则 转 (4), 否则 将 N 加 上 1, 用 
ve 代 wi 转 回 (2)。 
(4 追寻 最 短路 , 从 终点 n 开始 道 向 逐次 求 最 短路 经 过 的 点 (可 能 不 止 一 条 最 短路 ), 权 和 为 P(n)。 
3 算法 拓展 与 MATLAB 程序 实现 
求 n 个 顶点 的 连通 图 G 上 任意 两 点 间 的 最 短路 长 的 算法 ( 它 是 下 面 要 讲 的 Floyd(1962) 算 法 的 改 
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进 ,但 它 只 适用 于 非 负 权 的 最 短路 问题 )。 设 /为 顶点 i 到 j 的 距离 , 令 图 6 的 权 和 矩阵 DO = (dw) ,为 
lj, 当 (ij)eE 
d(0 =40, 当 i= 7 时 
,其 他 情况 


算法 基本 步骤 为 5: 

(输入 权 和 矩阵 Do. 

(2) 计 算 忆 = (dy)ws， 其 中 dy =min[d? +dy] 表 示 从 点 i 到 点 j 的 距离 , 或 者 经 过 某 一 个 中 间 点 到 达 j 
点 的 最 短路 。 

(3) 若 D= DO0 , 则 4 就 是 i 到 j 的 最 短路 长 ,不 然 , 取 P(0. =D, 转 (2)。 

其 中 第 (2) 步 最 多 执行 max1=floor(log(2*n-3)/log(2)) 次 。 其 原理 是 :第 (2) 步 得 到 的 矩阵 的 元 素 dj 表示 
从 i 点 到 j 点 的 距离 ,或 i 点 能 经 过 另 一 点 k 到 达 j 点 的 距离 中 的 最 小 值 。 因此 每 执行 一 次 步骤 (2) 所 得 
到 的 矩阵 D 相当 于 男 一 个 图 的 矩阵 ,这 个 图 的 构造 是 : 当 dd; 的 值 小 于 原 和 矩阵 的 值 时 ,将 原 图 上 将 雯 点 
之 间 加 上 一 条 权 为 dj 的 边 , 若 原来 有 边 , 则 改变 它 的 权 为 必 . , | 

设 ij 点 间 有 一 条 最 短路 ,这 最 短路 有 m 个 顶点 ,不 妨 记 这 些 顶 点 为 i it1yi+2,…jitm-1=j, 经 过 一 
次 (2) 的 步骤 时 , 求 得 (i, i+2), (i+2,i+4),…… 的 长 度 ,再 经 过 一 次 (2) 的 步骤 时 ,得 到 (i, i+4), (i+4,i+8),…… 的 
长 度 ，…… . 因为 当 2% >= m-1 >= 2^(k-1) + 工时 , 最 多 执行 有 次 可 以 求 得 点 1 到 点 m 之 间 的 最 短 距 ， 
所 以 重复 步骤 (2) 最 多 floor(oe Om oa RED] 的 最 短 距离 。 

编制 MatLab 程序 如 下 : 

function D=floyd(D0) %D0 是 原 关 联 和 矩阵 ，D 是 表示 最 短 距 离 的 矩阵 ,D(ij) 等 于 i 到 j 之 间 的 最 短路 
的 距离 : 

[n,nl=size(DO);D=zeros(n); % 创 建 和 矩阵 D, 虽 不 是 必要 的 ,但 这 样 程序 运行 得 较 快 

max]1=floor(log(2*n—3)/log(2)); 

for k=1:max1 

for r=1:n; D1=DO(r,:);  D(r,:)=min(repmat(D1",1,n)+D0O); end 

DO=D; % 此 前 可 加 上 判断 两 矩阵 是 否 相 等 的 语句 来 减少 循环 次 数 - 

end 

注 : 实际 上 步骤 (2) 执 行 的 次 数 常 常 不 到 maxl, 只 要 在 步骤 (2) 得 到 的 矩阵 D 与 DO 相同 即 可 停止 计 
算 , 因为 再 执行 (2) 仍 然 得 到 相同 的 矩阵 D. 这 时 , 可 在 程序 中 D0=D 之 前 加 上 判断 两 矩阵 是 否 相 等 的 
语句 。 但 因 和 矩阵 中 可 能 有 inf, 不 能 直接 进行 比较 。 具 体 语句 如 下 : 
DT=D; % 因 为 不 能 改变 D, 故 要 用 其 副本 DT 代替 D 
DT(ind(DT==inb)=1; Dodind(DO==inp)=1;DO=DO-DT; 
if sum(DO(:))==0 return; end 

以 上 程序 的 缺点 是 没有 标 出 最 短路 的 路 径 。Floyd 的 原 算 法 需要 计算 n 个 循环 ,但 是 可 以 确定 最 短 
路 的 路 径 , 且 适用 于 有 负 权 的 最 短路 问题 ,但 图 中 不 能 有 总 权 为 负 的 圈 。 算法 如 下 : 

(1) 输入 权 矩 阵 D(0) =(do )=D. 在 do 边 上 做 上 路 径 标记 (ij)， 

(2) 计算 DO = (dd 的) k=1,2,…,n, 其 中 d=min[dy",dx"+dy"] 表 示 在 第 步 时 从 点 i 到 点 j 的 某 路 
径 距 离 等 于 原 路 径 的 距离 和 经 过 中 间 点 k 到 达 j 点 的 距离 两 者 的 较 短 者 。 对 每 个 dd 做 上 路 径 的 标记 ， 
当 d = db 时 标记 不 变 , 不 然 将 dy 的 标记 为 qt 与 必 的 标记 合并 ， 去 掉 一 个 连接 处 重复 的 而 成 。 

(3) Dl 中 元 素 d0 就 是 i 到] 的 最 短路 长 ， 妹 应 的 标记 就 是 最 短路 径 。 

可 编 一 个 简单 Matlab 的 示意 程序 为 : 

function [C,D]=floydn(DO) 

[n,n]=size(DO); C=cell(n); D1=D0; D1(find(D1==inf))==0; 

M=2* max(D1(:)); [HIUJ]=find(D1); 

for ii=1:s C{HKiD,IGiiDj=fstrcat(int2str(IHiiD), 一 int2str(GJ(iD)))] end 

for k=1:n ;D=min(DO,repmat(DO(:,k),1,n)+repmat(DO(k,:),n, 1)); 

D1=D; DO(find(DO==inf))=M; D1(find(D1==inf))=M; DO=DO-D!1; 

[ILUI=find(DO>0); s=length(D); 
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for ii=1l:s_ C{1IGii,IU(Gii)j=strcat(C{Iii),kj,C{k,IJGii(2:engd)); end 
DO=D; end 
4 实例 分 析 


实例 1 使 用 最 短路 算法 求解 图 1 点 1 至 点 6 的 最 短路 。 
分 析 : 为 记号 简 清 起 见 , 省 略 了 步骤 数 N. 
(1) P(1)=0, 
(2) Ti(2) = P(1) + 7 = 7，Ti(3) = P(1)+9=9, 
(3) P(2) = min(T(2), Ti(3) = 7, TD = 1 
(2) Ti(3) = P(2) + 1 = 8，Ti(4) = P(2) + 7 = 14 
(3) P(3) = min(T(3), Ti(4)) = 8, T2(3)=2 
(2) Ti(4) = P(3) + 6 = 14 (也 等 于 旧 的 T,(4))， 
Ti(S) = Ti(3) + 6 = 14 
(3) P(4) = min(Ti(4), Ti(5)) = 14, P(S) = min(T;(4)，Ti(3)) = 
14, Tx(4) = (2, 3)，Tz(5) = 3， 
(2) Ti(6) = min(P(5)+6, P(4)+7) = 20， 
(3) P(6) = 20, Tz(6) = 5, 
(4) 最 短路 的 逆向 为 (6, 5, 3, 2, 1), 因此 最 短路 为 (1, 2, 3, 5, 6). 权 和 为 P(6)=20. 
实际 上 可 以 列表 1 简明 扼要 地 得 到 结果 : 


— 


表 1 中 括号 中 的 数 代 表 顶 点 ,第 
二 列 是 逐次 得 到 的 具有 P 标记 的 顶 
点 ,因此 第 二 列 第 一 行 是 起 点 (1), 第 
闵行 第 +1(>2) 列 括号 前 的 值 是 顶点 
(k) 的 TI((k), 和 NV) 值 ,括号 中 是 顶点 (k) 的 
Ts((k),N) 值 ,INF 表示 无 穷 大 (Infinity)， 
写 有 PP 的 表示 该 列 的 项 点 已 具有 P 
标记 , 不 再 参加 重新 标记 (在 实际 列表 
时 ,第 一 列 可 以 不 要 , 并 且 有 PP 的 单元 也 可 以 空 着 不 写 , 表 1 中 这 样 写 是 为 了 表达 得 更 清楚 ) 。 

实例 2 已 知 某 地 区 的 交通 网 络 是 一 个 连通 图 , 其 中 交点 表示 居民 小 区 , 边 表示 公路 , 边 长 为 权 ， 
权 和 矩阵 为 也 。 

0 30 inf inf inf inf inf 6 18 4 


inf inf 60 30 0 inf inf 15 0 
inf inf 25 18 inf 0 15 2 30 1 
inf 15 inf inf inf 15 0 图 2 


问 区 中 心 医院 应 建 在 那个 小 区 ,可 使 离 医院 最 远 的 小 区 居民 就 诊 时 所 走 的 路 程 最 近 ? 

分 析 :(1) 求 出 各 小 区 记 到 小 区 上 广 之 间 的 最 短 距离 的 和 矩阵 D. 当 边 是 无 向 时 ,D 是 对 称 的 ,小 区 立 到 
小 区 方 之 间 的 来 回 最 短 距 离 为 D2=D+D'=2*D; 当 边 是 有 向 时 , 小 区 立 到 小 区 方 之 向 的 来 回 最 短 距 离 为 
对 称 矩 阵 D2=D+D' ,其 中 D' 为 D 的 转 置 。 

(2) 行 向 量 far=max(D2) 的 第 广 个 元 素 表示 若 医 院 建 在 广 点 ,最 远 的 小 区 到 医院 的 来 何 距 离 。 

(3) [ m, 记 =min(far) 就 得 出 医院 应 建 在 认 点 ,最 远 的 小 区 到 医院 的 来 回 距 离 为 m. 

主 程序 为 : 

DO=inf*ones(7); 

DO(1,L1,2)D={0,30]; D0(2,[1,2,3,6])=[30.0,20,15]; 

DO(3,{2,3,4,5,6])=[20,0,20,60,251; 
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DO(4,[3,4,5,6])={20,0,30, 18]; D0(5,[3,4,5J)=[60,30,0]; 
D0(6,[2,3,4,6,7])=[1S,25,18,0,15]; 

D0(7,[6,7])=[15,0]; 

D=floyd(DO); D2=D+D'; far=max(D2); [m,ii]=min(far) 
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MATLAB Program of the Shortest Path Problem of Graph Theory 
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Abstract: This article puts forward the shortest path problem of Graph Theory, and sets forth the optimal algorithm——Di~ 
jkstra algorithm and its theory. Furthermore it describes and develops the model’s algorithm. Finally the article presents the MAT- 
LAB program of the shortest path problem. 

Key words: shortest path problem; dijkstra algorithm:; MATLAB program; shortest path length 


RR 


《上 接 第 20 页 ) a 
2 jc T(c+r+k)/(s + 

TR ON 

2 je T(c+r+p)/se” 


tk 


(9) 当 0<x<l1 时 ,有 P(R(1) <xlf,…,t,)= 一 [te dk 


两 边 对 * 求 导 ,可 得 R(1) 的 后 验 密度 为 


XM (Inx) "+Inx) 


hr (xh,...t,) = 下 ,0<x<l 


erat+r1 > 
k 


k=0 
在 试验 M ,,M ,,M ， 二 : R(t) 的 Bayes 估计 完全 类 似 地 可 求 得 。 


= 天 ctrtk 
-1D) IT(ctrt+k)/s 


参考 文献 : 

[1] 张 尧 庭 , 陈 汉 锋 . 贝 叶 斯 统计 推断 [M]. 北京 :科学 出 版 社 ,1991. 

[2] 黄 江 平 .定时 与 定数 截 尾 试 验 参 模 的 建立 与 参数 估计 []. 纯 粹 数学 与 应 用 数学 ,2004. 

[3] 张 丕 一 .定数 截 屁 寿 命 试验 下 指数 分 布 平 均 寿命 的 置信 限 的 优良 性 [. 青 岛 大 学 学 报 ,1999. 
[4] 丁 元 奖 ,. 指 数 分 布 寿命 试验 的 Bayes 可 靠 性 分 析 [J]. 宁 波 大 学 学 报 ,1998. 


Bayes Estimation of Exponent Distribution Reliability under Censoring Testing 
FANG Lian-di, GUAN Shi-jun | 
(School of Mathematics and Computer Science, Anhui Normal University, Wuhu 241000,China) 
Abstract: Under the exponential distribution fixed number or fixed time censoring, this paper discusses Bayes estimation of 
the R(t),when the prior distribution is different. 
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